Anti-phase synchronization of phase-reduced oscillators using open-loop control 
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In this letter, we present an elegant method to build and maintain an anti-phase configuration 
of two nonlinear oscillators with different natural frequencies and dynamics described by the sinu- 
soidal phase-reduced model. The anti-phase synchronization is achieved using a common input that 
couples the oscillators and consists of a sequence of square pulses of appropriate amplitude and 
duration. This example provides a proof of principle that open-loop control can be used to create 
desired synchronization patterns for nonlinear oscillators, when feedback is expensive or impossible 
to obtain. 
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Natural and artificial systems that consist of collec- 
tions of isolated or interacting nonlinear oscillators are 
reaching complexity levels that are beyond human un- 
derstanding. The normal operation of these complex 
systems often requires the formation of certain synchro- 
nization structures. The synchronization of oscillating 
systems is a fundamental and extensively studied phe- 
nomenon in natural sciences and engineering [l|, Q • Ex- 
amples include circadian rhythms [3j|, neural circuitry in 
the brain Q, pacemaker cells of the heart Q, insulin- 
secreting cells of the pancreas @, chemical oscillations 
Q, semiconductor lasers in physics Q, and vibrating 
systems in mechanical engineering Q. The above sys- 
tems are often extremely large in scale, which poses seri- 
ous theoretical and computational challenges to efforts to 
model, guide, or optimize them. Deriving control signals 
that can drive complex systems to desired synchroniza- 
tion configurations is of utmost theoretical and practical 
importance [lfj| . A premier example comes from the area 
of neuroscience, where devising low-power external stim- 
uli that synchronize or desynchronize a network of cou- 
pled or uncoupled neurons is imperative for wide-ranging 
applications, from neurological treatment of Parkinson's 
disease and epilepsy |ll| - fl3j to the design of neurocom- 
puters 0, \M- 

Mathematical devices are required for describing the 
complex dynamics of oscillating systems in a manner 
that is both tractable and flexible in design. A promis- 
ing approach to construct simplified yet accurate mod- 
els that capture the essential overall system properties is 
through the use of phase model reduction. Underlying 
this method is the fact that an oscillating system with a 
stable periodic orbit, described by a system of first or- 
der ordinary differential equations, can be reduced, under 
certain circumstances, to a dynamic equation in a single 
variable, whichjepresents the phase 9 of the system on 
its limit cycle 
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where u(t) is some external input, for example a current, 
f(9) gives the baseline dynamics, which is present even 
when u = 0, and Z(ff) describes the phase sensitivity 
to the stimulus and is called the phase response curve 
(PRC). For example, for the prototype mathematical 
neuron model, the famous Hodgkin-Huxley model which 
describes the propagation of action potentials along the 
giant squid axon and is a complex system of four highly 
nonlinear differential equations the corresponding 
phase-reduced model can be computed [l7j using widely 
available software |18j |. 

Phase models have been very effectively used in theo- 
retical, numerical, and more recently, experimental stud- 
ies to analyze the collective behavior of networks of oscil- 
lators 191-422 [ . Specifically, these models have been used 
to investigate either the patterns of synchrony that result 
from the type and architecture of coupling [n], 23 T -26 1 . 
or the response of the system to external stimuli 1J, |27| ■ 
Motivated by studies in control and systems theory which 
show that feedback is essential for investigating properties 
of complex dynamical systems such as self-organization 
and stability [28j |. various feedback approaches, based on 
phase models, have recently been developed for the ef- 
ficient control of synchronization patterns in oscillator 
assemblies [l(| HH, H3 ■ Although these synchronization 
engineering methods are effective for the synthesis of sub- 
tle dynamical patterns such as itinerant cluster dynam- 
ics, in many emerging applications involving the control 
of large-scale complex systems, state feedback maybe dif- 
ficult, impossible, or expensive to obtain, or the types of 
feedback laws that can be used are severely restricted due 
to the complexity of system dynamics. As a result, the 
development of open-loop controls for the design of phase 
patterns for oscillator ensembles is compelling. 

Recently, geometric control theory |31| has been em- 
ployed to study controllability of a network of oscillators 
with different natural oscillation frequencies 32]. Con- 
trollability guarantees the existence of open-loop inputs 
that can drive a system of oscillators to any desired syn- 
chronization pattern. The idea to use an open-loop con- 
trol to achieve the target synchronization structure is still 
useful even when the system is not fully controllable but 
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the final state is reachable 

In this letter, we demonstrate the potential of open- 
loop control to engineer patterns of synchronization in 
collections of nonlinear oscillators, using a simple yet im- 
portant system. Specifically, we show how an input se- 
quence composed of surprisingly simple square pulses can 
be used to build and maintain a 7r phase difference be- 
tween two nonlinear oscillators with dynamics described 
by a sinusoidal phase-reduced model, with / constant and 
Z(Q) = sin0 in (JXJ) . This model is an approximation to 
the FitzHugh-Nagumo model, a two-dimensional simpli- 
fication of the Ho dgkin- Huxley neuron model, near the 
Hopf bifurcation point [3J]. If Parkinson's disease is con- 
sidered as the synchronized response of an assembly of 
neurons, described by phase models, then it is desirable 
to bring these oscillators out of phase with an open-loop 
signal. This observation highlights the importance of es- 
tablishing an anti-phase configuration for a system of two 
phase-reduced oscillators. 

In the following we consider a pair of oscillators de- 
scribed by the sinusoidal phase-reduced model 



ix=bj\ + usmfti, 
9 2 = cj 2 + usin0 2 , 
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where u; 2 > ^1 but close to each other. The common in- 
put u(t) couples the two oscillators and we suppose that 
it is bounded by |u| < 0J1. We assume that a well-defined 
initial state can be established by a classical phase re- 
setting method, where a large pulse is applied to the 
oscillators and brings them to the same initial phase, in- 
dependent of the phases they had before @. The goal 
is, starting from 0i(O) = 62(0) = 0, to build a ir phase 
difference between the two oscillators and then maintain 
this anti-phase configuration. We describe how this can 
be done efficiently using simple square pulses. The ap- 
proach is different from that in where the phase 
difference is built asymptotically in time. 

First we show how this anti-phase configuration can 
be produced by applying a constant signal u = —M of 
appropriate magnitude for a duration Tm- If we require 
at time t = Tm that 

0i (T M ) = (2n - 1)tt, 2 (T m ) = 2nn 

for some positive integer n, then, using ^ and ([3]), we 
obtain 
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or, after integrating, 
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Let L(M) and R(M) denote the left and right hand 
sides of ((H), respectively, as functions of the amplitude M. 
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then L(0) < R(0), while it is always L(lo\) > R(u)i). 
Thus, for lo\ sufficiently close to uj 2 , the transcendental 
equation ((4} has a solution M S (0,wi). The duration 
Tm of the constant signal u = — M is given by 



T M = R{M) 
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Note that in the absence of any input, i.e. when u — 
0, a time Tq = Tr/(u>2 ~ ^1) is necessary to build a ir 
difference between the oscillators. Also, observe that for 
given frequencies w\ and W2, © can be used to calculate 
the allowed values of n. As n increases, we expect that a 
signal with lower amplitude M and larger duration Tm 
is required to create the 7r phase difference. 

Once the anti-phase configuration has been set up, we 
can maintain it by appropriately choosing the control u. 
The idea is that when 2 ~ $1 — ^ , then sin0! and sin0 2 
have always opposite signs, so we can use as common 
input a square wave with appropriate amplitude m and 
period T m 



u(t) = (-l)V T M + 
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k = 1,2,..., to accelerate the slow oscillator and decel- 
erate the fast one. 

The control amplitude m is chosen to assure that phase 
difference is preserved at the switching instants of the 
square wave 
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This implies that the modified half-periods of the two 
oscillators for t > Tm , after the application of the square 
wave, are the same and equal to the half-period T m /2 of 
this wave. If, without loss of generality, we use the first 
half-period Tm < t < Tm + T m /2 after the onset of the 
square wave, we obtain the condition 
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or, after the change of variables 4>i = $1 — (2n — l)7r, (j>2 — 
2 - 2wr, 



ui% + m sin (pi J ui2 — rn sin <P2 



d<t>2 



(7) 



Observe that different integration intervals in (|6]), due to 
the already built 7r difference, led to a different control 
sign in the denominator of the integrands in (J7J, increas- 
ing thus the period of the fast oscillator and decreasing 
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FIG. 1. (Color online) For cji = 1.9-7T rad/sec and CJ2 = 2.17T rad/sec we plot the input function (a, c, e) and the phase difference 
02 — 6i (b, d, f), for n — 1,3,5. Observe that for larger n, a smaller amplitude is necessary to build the n phase difference. 
After this difference has been built, a simple square wave is sufficient to maintain it, the same for all n. The small periodic 
deviation from ir is expected since the method is designed to preserve the anti-phase configuration at certain instants separated 
by half period of the square wave. 



that of the slower. Integrating ([7]) yields 
2 




(8) 



If we use L(m) and R(m) to denote the left and right 
hand sides of (|8]), as functions of the amplitude m, then 
L(0) = 7r/wi,i?(0) = tt/w 2 and L{uj x ) -> 2/u) 1 ,R(u 1 ) > 
since uj 2 > u>i, it is L(0) > R(0). The 
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requirement L(oj\) < R(uii) is satisfied when 2/wi < 



7r/Va 



i\, which leads to the condition 
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When © holds, the transcendental equation ([5]) has a 
solution m € (0,cji). The modified period is 



Tn 



2L(m) 



: tan 
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For = 1.97T rad/sec and W2 = 2.17T rad/sec, con- 
dition ((9]) is satisfied and from ([5]) we find n < 5. In 
Fig. Q]we plot the input sequence (left column) and the 
resulting phase difference (right column), for n = 1,3,5. 
Observe that for increasing n the necessary amplitude 



M to build the desired phase difference is smaller but 
the corresponding time Tm is larger. The maintenance 
signal is the same in all cases. Note that after the estab- 
lishment of the phase difference there is a small periodic 
deviation from 7r, but this is expected since the method 
is designed to maintain the anti-phase configuration at 
certain time instants separated by At = T m /2. 

We have demonstrated that open-loop control can be 
used to design patterns of synchrony for nonlinear oscil- 
lators, by presenting a simple and elegant method which 
can achieve and maintain a tt phase difference between 
two phase-reduced oscillators using surprisingly simple 
square pulses. This approach can find applications across 
various disciplines in natural sciences where synchroniza- 
tion plays an important role, to situations where feedback 
is expensive or impossible to obtain. 

This work was supported by the Air Force Office of 
Scientific Research under grant #FA9550-10-1-0146. 
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